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The applications of elastic and quasielastic light scattering techniques to polymers in dilute solution, ther- 
motropic liquid crystals and biological membrane vesicles are presented. It is focused on how we extract specific 
structural features or dynamic processes of these condensed medium samples through the light scattering 
methods. The paper reviews the light scattering studies that were carried out at the author's laboratory in 
Wisconsin from 1973 until 1980. 
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1. Introduction 

With the advent of the dynamic light scattering method 
and the attendant detection techniques, [1-5], we are now 
able to expand the application of light scattering methods 
to a variety of diverse condensed medium systems hitherto 
unexplored. In this paper, I will outline the kinds of light 
scattering methods that we use in my laboratory at Wiscon- 
sin to probe the structure and dynamics of polymers, 
nematic liquid crystals, and biological membranes. Four 
scattering techniques shall first be described with ap- 
propriate examples in each case, and then I will move on to 
the studies of (a) amphoteric latex system, (b) photo-receptor 
disk membrane vesicles, (c) binary nematic solution, and (d) 
intrachain dynamics of random coil polymers, in order to 
bring home the power and limitation of various light scat- 
tering techniques. At the outset I must emphasize that this 
paper is not intended as an exhaustive review of the state-of- 
the art of light scattering methods, but rather as a report of 
how one academic research laboratory uses these methods 
to explore diverse problems of interest. 

2. Elastic Light Scattering 

Extensive treatises have been written on the subject such 
that I need not dwell on the historic overview or lengthy ex- 
plication of the technique. Instead, I will focus on a rather 



narrowly defined problem of how to determine the linear 
dimension of scattering particles with substantial symmetry 
in shape, particularly when the linear dimension R is com- 
parable to the incident wavelength X. The problem of this 
kind arises for example in trying to determine the large 
radius of gyration [9] of T even bacteriophage DNA whose 
molecular weights are on the order of 10 8 daltons. If one 
tries to effect the customary procedure of the Zimm plot [10] 
to extract the radius of gyration, the scattering intensity 
profile must be obtained at such small angles that one en- 
counters substantial technical problems [9]. By small angle 
scattering, I mean that one must obtain the data in the 
Guinier region, [11] i.e., QR < 1 where Q is the momentum 
transfer. We propose a new scheme to determine the linear 
dimension R when the scattering particles are so large that 
QR < 1 is difficult to attain experimentally. The method is 
to focus on the structure of scattering profile at higher 
angles where the Bragg condition is fulfilled, i.e., QR = 
m/ir where m is an integral multiple constant. To be more 
specific, I quote some simple examples such as solid sphere, 
hollow sphere, spheroidal shape of either oblate or prolate 
axial ratio, and the corresponding shell structures. The 
isotropic parts of the particle form factor in the limit of 
Bayleigh-Gans-Debye scattering [12] of some of these are 
given below. 



a. Solid sphere 



1 Figures in brackets indicate literature references at the end of this paper. 



P(0) = \jJLx) + j 2 (x)Y 



(i) 
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where x = QR, R is the radius and j k (x) is the &th order 
spherical Bessel function. 

6. Hollow sphere [13, 14] 

P(9) = [ — y . (sin* — sinxl — *cosx + z/cos:r/)J 2 (2) 

where x = QR, I = aJR, R and a are respectively the outer 
and inner radii. 



c. Ellipsoids of revolution [15] 



m = 



9tt 



Jhi(u) 



cos(3d(3 



(3) 



where u s Qa(cos 2 j3 + -=■ sin 2 /3) H , a and 6 are semi-major 
a 

and semi-minor axes, (8 is the angle between the semi-major 

axis and bisectrix, and J 3/2 is the 3/2th order Bessel 

function. 



d. Circular cylinder 
'2v 



wsir-ssr-^c-ffl-^ai'-^w 



where a and / are the radius and length of the cylinder, u = 
Qa, v = Ql/2, and /3 is the angle between the cylinder axis 
and bisectrix. 



e. Spheroidal shells 



m= r' sin'fWl -q* l dt 
K r>(l - qt 2 ) 

_ ,' sin 2 [y(-t^f-)H] 



/(-¥=£-) 



dt 



(5) 



(6) 



for oblate shell where x = Qb, y = Qa, q = 1 - (-ff, and 
6 and a are respectively the semi-major and semi-minor axes 
such that < q < 1. 



m= r'sin 2 t*yi +pt 2 J d 
K ^(1 + pt 2 ) 



sin 2 \y{- 



1 -*• •"' \u 



.)Vi] 



f(-¥£-) 



(7) 



Jf (8) 



for prolate shell where x and y are the same as above, a and 
b are respectively the semi-major and semi-minor axes, and 
p s (a/6) 2 - 1 such that < /? < oo. 

As an illustration of the proposed method, we take the 
simplest structure, i.e., sphere, and show how the Bragg 
condition is extracted from which the radius of sphere is 
deduced. The isotropic particle form factor of sphere in the 



Rayleigh-Gans-Debye limit is given by the square of the sum 
of the zeroth and second order spherical Bessel function as 
in eq (1), and its structure is a monotonically decreasing 
function from unity at zero scattering angle to x = tana; 
where it gives the first minimum. This is illustrated in 
figure 1 where a semi-logarithmic plot of the particle form 
factor F\6) versus sine of one-half of the scattering angle 6 is 
given for spherical particles having the radii of 100, 110, 
120, and 130 nm suspended in water and scattered by inci- 
dent blue light of 436 nm in wavelength. It is clear from the 
plots that under these conditions the radius must exceed 
120 nm for P(d) to give the first minimum; P(d) is fairly 
structureless for the radius less than 120 nm. If, on the other 
hand, we plot (QR)*P{6) against sin(0/2), the graphs as shown 
in figure 2 result. This arises because the damping profile of 
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FIGURE 1. The isotropic part of the particle form factor P(x) of Ray- 
leigh-Gans-Debye scattering from solid, isotropic spheres against sin (0/2) 
for the radius of 100, 110, 120 and 130nm, with the incident wavelength of 
436 nm (blue) in a medium with the refractive index n of 1.33 (water). Note 
that x= QR hence P(x) = P(e). 
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Figure 2. x*P(x) versus sin (0/2) of figure 1. In place of R = 130 nm, 
plot for 200 nm is shown. 
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P[6) in eq (1) is compensated by the factor (QR)*. Thus, the 
monotonic descrease of P(0) for small particles is rendered 
to have a maximum by multiplying by (QR) 1 . Algebraically 
this is easily seen from eq (1); 



** P(d) = 9 (- 



— cos*) 2 



(9) 



Airn . 



where x = QR = -r — sin{dl2)»R and the extrema posi- 

Ao 

tions of P(6) appear at 



tan* = 



_ x 



1-x 2 



for maxima 



and 



tan* = * for minima 



The plots in figure 2 are drawn for the common experimen- 
tal scattering conditions in mind, namely the scattering 
angle spans from 30° to 150°. Experimentally the ordinate 
scale is immaterial to the extent that one is interested in 
determining the extrema positions to extract the particle 
radius. For completeness sake, we show in figures 3 and 4 
the analogous set of plots for spherical particles with longer 
radii as though they could be analyzed by the Rayleigh- 
Gans-Debye approximation, whereas its applicability is 
limited [8] in this size range, and P(0) itself contains suffi- 
cient structure that x*P(d) plot affords at best sharper max- 
ima and smoother minima. I shall later return to the limit of 
applicability of Rayleigh-Gans-Debye scattering relative to 
the particle size. We display in figure 5 why the factor ** is 
chosen as the multiplication factor to P{6) in order to 
moderate its damping at higher angles. For a given radius 
of a particle, we show three profiles whereby it is made clear 
that a factor x 2 does not reduce the damping sufficiently, 
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Figure 3. Similar to figure 1 with R = 200, 300, 400 and 500 nm. 



whereas a factor x 6 overcompensates the damping. The op- 
timum is clearly the x* factor. 

Turning to the experimental verification of the proposed 
method, [16], we show in figure 6 a plot of I vv (6)»siti*(dl2) 
versus sin(0/2) for a polystyrene latex standard (Poly- 
sciences, Lot 2-1435), whose radius is given as 87 ±4 nm 
determined by transmission electron microscope, sus- 
pended in distilled water. Here, the subscripts vv of the scat- 
tered intensity stand for the polarized scattering, namely 
the polarizer and analyzer are both oriented vertically 
relative to the scattering plane. We should note that the or- 
dinate scale is arbitrary and x*P(d) is proportional to 
I„(6)»sm*(dl2), since P{6) is the isotropic part of the particle 
form factor. The data points were taken at 1 ° increments 
and the solid curve was drawn according to eq (9) with Xo = 
436, n = 1.333, R = 80 nm. The discrepancy of 10 percent 
between the radius of 87 nm by electron microscopy and 
that of 80 nm by this method may be ascribed to a number 
of artifacts arising from electron beam optics in the electron 
microscopy technique. In figure 7, we show a test of another 
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FIGURE 4. Similar to figure 2 with R = 200, 300, 400 and 500 nm. 
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FICURE5. x^i^fl) and xfP(6) versus x. 
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polystyrene latex standard (Polysciences, Lot 2380), having 
the radius of 139 nm again determined by electronic 
microscopy. Two curves are drawn for comparison. The first 
(solid curve) is that predicted by the Mie scattering function 
with R = 136 nm and the other (dashed curve) is that pre- 
dicted by eq (9) with the same R. In the Mie function fitting, 
we use the refractive index ratio m of the particle to 
medium as 1.21 [17]. In the inset, the scattering profiles at 
different concentrations of the latex particles are shown 
where the concentration range of the most (a) to the least (d) 
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Ficure 6. Calibration run of the light scattering method with Polysciences 
latex standard with R = 87±4 nm (mean ± s.d.). The scattering intensity 
here is observed with the incident and scattered beams vertically polarized 
relative to the scattering plane. The ordinate is scaled in arbitrary units, 
and the solid curve represents the Rayleigh-Gans-Debye scattering func- 
tion for an isotropic solid sphere with 80 nm raduis. 




FiCURE 7. Calibration run with Polyscience latex standard with R = 135± 
1 nm. Experimental conditions are identical with those in Figure 6. Dashed 
curve is drawn with the Rayleigh scattering function with R = 136 nm and 
solid curve with the Mie function with R = 136 nm. The inset shows the 
concentration dependence of scattering profile whereby a represents the 
most concentrated and d the least concentrated. The number density of 
latex particles in a is about one order of magnitude larger than in d. 



concentrated is about one order of magnitude. It shows well 
that the scattering profiles are weakly dependent on the 
concentration, and we could have deduced the radius of 130 
nm from the most concentrated case (a) shown in the inset. 
The least concentrated case (d) is expanded in the plot of 
figure 7. Ignoring the matching with the entire scattering 
profile but focusing on the maximum and minimum posi- 
tions inferred from eq (9) for the Rayleigh-Gans-Debye scat- 
tering, we would have deduced the radius of 144 nm. Thus 
for this size particle, either scattering function would have 
sufficed if the radius determination within 4 percent is 
acceptable. It should however be noted that the discrepancy 
of the radius determinations by the two scattering functions 
exceeds the precision limit of extrema position determina- 
tion in this size range. 

Having thus established that our proposed method works 
for the radius determination of a spherical particle of R ^ 
140 nm suspended in aqueous media, we now turn to the ap- 
plicability limit of Rayleigh-Gans-Debye scattering and 
when one must use the Mie scattering function. In order to 
stipulate the applicability limit, we compare the extrema 
positions predicted from eq (9) and those computed from 
the Mie functions. The latter computation was performed 
on a Harris/7 computer. The comparison is provided in 
figures 8(A) and 8(B) where the reduced size parameter a, 
defined as 2ttR/K, is given in terms of the extrema positions 
in the scattering angle of X*F\6) profiles for the polarized 
and unpolarized scatterings respectively. When m = 
1.0001, we recover the predictions of the Rayleigh-Gans- 
Debye scattering as given by eq (9); this particular value of 
m is not significant as long as m is very close but still larger 
than unity for the Mie function to be evaluated. These 
figures allow us to understand why the radius deduced from 
the Rayleigh function always overestimates. At any given 
extreme position, be it a maximum or minimum, a values 
for m > 1 invariably lie below that for m = 1 although the 
relative error committed by assuming m « 1 is not a mono- 
tonically increasing function of m for any size particles. 
This is particularly true with the maximum positions which 
show oscillation with respect to m. Hence, the relative error 
should also oscillate. In fact, we can evaluate the error. This 
is illustrated in figure 9(A) and 9(B) where the percent error 
of analyzing the scatterng profile extrema according to the 
Rayleigh-Gans-Debye scattering is plotted against the 
reduced size parameter a at different values of m. The 
polarized and unpolarized scattering cases are shown in (A) 
and (B), respectively, In either case, the error can be equal 
or less than 10 percent if m < 1.15 for a < 7. The oscilla- 
tion in the error estimate becomes progressively larger in 
amplitude as the m value increases. For the values of m 
close to unity, the error with respect to the particle size is 
sensibly constant though with a slight increase with a in 
each m. Hence, one should not approximate with impunity 



574 



the scattering of small particles by the Rayleigh scattering 
function regardless of the refractive index ratio m. 

In closing this subject, let me emphasize that the spher- 
ical particle size analysis can be effected by determinations 
of the extrema positions of sin 4 (0/2)/(0) profiles and these 
are provided by analytical solutions of 



x = tanx 



(minima) 



1 - x 2 



= tan* (maxima) 
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FIGURE 8A. Ratio of isotropic solid sphere's circumference to wavelength 
in scattering medium, a= 2nR/X, versus sin (0/2) at different extrema 
positions of the Mie function for x*F\8). Each set is drawn at different 
refractive index ratios m. The maxima and minima are distinguished by 
solid and dashed curves respectively. This is for the Mie function at the 
vertical/vertical optical configuration. 
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FIGURE 9A. Percent error of analyzing the profile extrema positions 
according to Rayleigh-Gans-Debye scattering as a function of the reduced 
size parameter a at different refractive index ratios m, for the polarized 
scattering. 
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FlCURE 8B. Analogous to (8a) for the case of Mie function at unpolar- 
ized/unpolarized configuration. 
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FIGURE 9B.Analogous to (9a) for the unpolarized scattering. 



575 



in the case of Rayleigh-Gans-Debye scattering whereas they 
must be evaluated numerically in the case of Mie scattering. 
Because the method depends on the structure of scattering 
function at high Q, it is restricted to monodisperse systems 
and extrapolation to infinite dilution is not an essential step 
in the procedure provided the scattering suspension is 
dilute enough. 



3. Quasielastic Light Scattering 

Here, we restrict our discussion to the thermally induced, 
spontaneous concentration fluctuations whereby the trans- 
lational diffusion coefficient of the scattering particles at in- 
finite dilution is deduced. The homodyne power spectrum 
[1] S(Q,v) of the Doppler broadened scattering from a 
monodisperse system of particles is 



S(Q,p) = A 



AVXA 



v 2 + (A^) 2 



+ B 



(10) 



where A is an optical constant which depends on the inten- 
sity factor of the spectrum, B is a constant, a measure of 



shot-noise level and the spectral half-width at half-height 
Api/i is related to the translational diffusion coefficient D by 



Ai»vi = DQ 2 Itc 



(11) 



Our instrument [18] is schematically depicted in a block 
diagram in figure 10. A typical power spectrum obtained 
from a Dow Polystyrene Latex standard suspension (45.4 
nm radius) is shown in figure 11 and the corresponding 
spectral halfwidth against Q 2 It is displayed in figure 12 
where the scattering angle spans the range 10°-100°. The 
diffusion coefficient deduced via eq (11) is (5.39 ±0.04) X 
10" 8 cm 2 /s which is in turn converted to the Stokes radius of 
45.5 ±0.4 nm. Progression of the S/N ratio of the observed 
power spectrum with the number of accumulations is shown 
in figure 13. 

4. Electrophoretic Light Scattering 

The technique was first developed by Ware and Flygare 
in 1971 [19] and subsequently by Uzgiris in 1972, [20] and it 
has since undergone substantial refinements [21-29]. The 
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F.CURE 10. Block diagram of Rayleigh spectrometer. 
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advantages that it offers over the conventional electrophor- 
esis methods have now been well-documented in the liter- 
ature [25, 27]. As with any new technique, this one also had 
to first be calibrated against those of more conventional 
methods with the use of a test system. Ware and Flygare 
chose bovine serum albumin (BSA) for the purpose because 
it was one of the best characterized globular proteins and 
commercially available in a relatively pure form. Others 
have subsequently chosen BSA for the same reason to 
calibrate their instruments [28-32]. 

The method is no more than another application of laser 
velocimetry. A monodisperse system of charged particles in 
dilute solution under the influence of an applied electric 
field would drift uniformly to the oppositely charged elec- 
trode. The diffusion equation governing this situation [4] in 



terms of the self part of the space-time autocorrelation func- 
tion G,(R,f) is 



-f- G M (R,t) + \ d .VG,(R,t) = Z)V 2 G,(R,0 (12) 
at 



with the initial condition 

GXR,0) = 5(R) (13) 

where V d is the uniform drift velocity, D the translational 
diffusion coefficient of particle and S(R) is the Dirac delta 
function. Upon taking the space Fourier transforms of the 
above, we have 

jrF.(Q,t) + iQ-Y d F,(Q,t) = - Q 2 DF,(Q,t) (14) 
o t 
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FIGURE 11. Homdyne power spectrum of a polystyrene latex standard at 2 



kHz bandwidth with Av 1/S = 299 Hz. 
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Figure 12. Av„ versus q 2 /n whose slope yields D = (5.39±0.04) X 
10~ u mVs and the Stokes raduis of 45.6 ± 0.4 nm. 



with 

F,(Q,0) = 1 (15) 

The solution of eq (14) with the initial condition, eq (15), is 

F,(Q,t) = e- Q2m e'°-V (16) 




FIGURE 13. Power spectra of polystryene latex standard taken at 60° scat- 
tering angle. A, B, and C are 250, 1000 and 4000 scan averages, respec- 
tively. 
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The power spectrum of heterodyne beating from a system of 
charged particles governed by eq (16) is a Doppler shifted 
Lorentzian whose shifted angular frequency Ao, is given by 



Aa>, = Q.V d 



(17) 



Our scattering geometry is represented in figure 14. In 
order to effect heterodyne beating, the scattering angle has 
to be fairly small (2°-8°) so as to make use of cell surface 
reflected light as the local oscillator component. Thus the 
shifted frequency 



Av, = QV d cos(0/2)/27r ~ Tmn(dV pp )l360d\ o (17) 



where n is the refractive index of scattering medium, n the 
electrophoretic mobility (the drift velocity per unit field 
strength), 6 the scattering angle, V pp the peak-to-peak ap- 
plied voltage, d the electrode spacing (1.84 mm in our case) 
and X„ the incident wavelength in vacuo. Hence, the shift 
frequency Av, should be proportional to the product, 6V PP , 
and the electrophoretic mobility /* can be deduced by deter- 
mining Av, measured at different scattering angles and V pp , 
Since there were some discrepancies in the mobility 
values of BSA reported by Ware and Flygare [19, 21] and 
Mohan et al. [28] by electrophoretic light scattering and 
those by Schlessinger [30], Alberty [31], and Longsworth and 



Jacobsen [32] with the moving boundary method, we have 
set out to examine the discrepanices [33]. The purpose was 
to test the accuracy attainable by this technique vis-a-vis 
that of a more conventional method. In so doing, we have 
not only found the mobility by this technique to be in com- 
plete accord with those of the moving boundary method but 
also established that the BSA mobility depends on the ionic 
strength of the suspending medium according to Henry's 
formulation [34, 35]. The latter finding is neither without 
parallel nor unexpected [36-38], but the ease of the experi- 
ment to confirm it points to the utility and power of the 
technique. At the same time, our experiment, which covers 
a wide range of ionic strength, points to certain limitations 
and the complementary nature of the technique with more 
conventional methods. 

The discrepancies referred to above were entirely due to 
the sample polydispersity and had little bearing on the 
veracity of the electrophoretic light scattering technique. 
This was confirmed by examining the two sets of samples. 
First is the so-called Fraction V of Armour (Lot A21505) 
which was used by Ware and Flygare as well as by all others 
with the moving boundary method, and the other is the BSA 
Monomer Standard of Pentex brand from Miles Labora- 
tories. Polyacrylamide gel electrophoresis patterns of the 
two are compared in figure 15. In figure 16, we display how 
well the shifted frequency depends linearly on the product 
(Q*V app ) for the Fraction V sample. Finally we compare the 
deduced mobilities, corrected to 20 °C in water, of the two 
samples relative to their ionic strength dependence in 
figure 17. In the same figure, we plot all other available 
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Armour Fraction V 




FIGURE 15. Polyacrylamide gel electrophoresis patterns of polydisperse 
BSA Fraction V of Armour and Pentex brand BSA monomer standard. 



BSA data in the literature at pH 8.7-9.4 after making ap- 
propriate temperature correction. Two solid curves are 
drawn according to Henry's equation 



* = i£fc«* R) 



(18) 



where Q E is the net charge of the electrokinetic unit, rj the 
viscosity of the medium, R the radius of the unit, x' 1 is the 
Debye screening length, and f{xR), which accounts for the 
ionic strength dependent charge screening, is given by 



J(xR)=(\+xRr\\ + 
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where the integral in the last term is the exponential in- 
tegral Ej[xR) [39]. With use of the Stokes radius of 36 A 
deduced from the translational diffusion coefficient of BSA 
by Baldwin et al. [40] and with the dependence of x~ l on 
ionic strength / at 20 °C as x~ l =3.045/ y/I(A), the two 
curves are drawn with 20 and 26 electronic charges per BSA 
molecule for Q E , respectively, for the lower (Fraction V) and 
upper (Pentex) sets. Since the data by the moving boundary 
method are all obtained with the Fraction V samples, it is 
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FIGURE 16. Electrophoresis Doppler shift frequency Av, versus 6V PP for 
Fraction V samples at 5 mM ionic strength. Three independently prepared 
samples are examined. Scattering angles of two preparations are distin- 
guished as (•) 2°, (3) 3°, (O) 4°, (6) 4° and third preparation is meas- 
ured at (□) 6°. 
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FIGURE 17. Electrophoretic mobility, corrected to 20 "C in water, versus 
ionic strength. Open circles with error bars are Fraction V sample, filled 
circles with error bars are Pentex sample. * is from the slope of figure 16. 
All others are from the literature; (□) Schlessinger, (O) Alberty, (O ar| d V) 
Longsworth and Jacobsen, (A) Ware and Flygare [20] and (A) Ware and 
Flygare22 and(O) Mohan et al. [29] 
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not surprising that they belong to the lower set and are con- 
sistent with ours. The results in figure 17 represent the first 
instance, to the best of our knowledge, of a systematic study 
of BSA mobility dependence on ionic strength at a given 
pH. It is clear that Henry's equation can well account for 
the ionic strength dependence of electrophoretic mobility. 



A model of dilute polymer solution is composed of opti- 
cally isotropic solvent and linear flexible chains, each of 
which is constituted of n+1 anisotropic elements with cylin- 
drical symmetry. The FDS spectral profile is formally given 
by 

I Vh (o>)=B(N/V) J_-. < £ t < (r)<>(0)> e-'^dr (20) 



5. Forward Depolarized Scattering (FDS) 

The technique was first proposed and utilized by Wada 
and coworkers in 1969-70 [41, 42] for the determination of 
the rotatory diffusion coefficient of a rod-like molecule, i.e., 
tobacco mosaic virus. A further refinement was reported by 
Schurr and Schmitz in 1973 with the same system [43] and 
an extension to calf-thymus DNA followed [44]. I will illus- 
trate the technique with use of an example, namely the 
intrachain dynamics of linear flexible macromolecules in 
dilute solution. A theoretical formulation of the spectral 
profile of a scattered optical field is provided by Ono and 
Okano. We have slightly generalized the scheme by con- 
cluding that the FDS spectrum is a superposition of uni- 
formly weighted multiple Lorentzian as long as there exists 
the normal coordinate transformation for the intramolecu- 
lar chain dynamic modes [45]. 



where B is an optical constant, (N/V) is the number density 
of polymer molecules within the scattering volume, the sub- 
scripts Fand h refer to the vertical and horizontal polariza- 
tion directions of the incident and scattered optical fields, 
respectively, and < > refer to the equilibrium ensemble 
average. If the normal coordinate transformation is possi- 
ble, the long wavelength modes in hydrodynamic regime is 



1/t* 



iytM-W/VHAaW & ^- Tkf 



(21) 



where Aa is the optical anisotropy of each element and t* is 
the kth normal mode relaxation time. In case of a monodis- 
perse system of dilute rigid rod molecules, the heterodyne 
power spectrum corresponding to eq (21) is 



S (0 v)-J (1/27rT) 
^'^-"V +(1/2x7)* 



(22) 
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FIGURE 18. Block diagram of the Rayleigh spectrometer in the FDS configuration and the filter train. Here, 1°, F and P stand for the incident, leakage 
and scattering intensities, respectively, and the subscripts V and H designate vertical and horizontal polarization directions. 
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where r is the relaxation time of rotatory diffusion of the 
rod. 

The optical train of our apparatus [46] is shown in figure 
18, and an example of S ± (0, v) for a dilute solution of 
poly(n-hexyl isocyanate), a rigid rod chain of 3500 A in 
length, is given in figure 19. 
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FIGURE 19 FDS power spectrum of PHIC (M.W. 3.2 X 20 s ) in n-hexane at 
C = 2.16 mg/mL with single Lorentzian fit ( — ) and Zimm spacing fit 
( — • — ). Upper and lower graphs are distinguished by linear and 
logarithmic frequency scales. 

The four light scattering techniques that we employ to 
study a variety of scattering systems are summarized above. 
Before leaving this section, I note parenthetically that our 
data acquisition and analysis system makes use of a mini- 
computer (PDP8/e) and a set of two microprocessor com- 
puters (Apple II) interacting with a Harris/7 system. This is 
shown schematically in figure 20 where only the part deal- 
ing with computer controlled automatic goniometer on 
SOFICA light scattering photometer is yet to be imple- 
mented while all others are now in operation. 



6. Systems 

Turning to the four systems that were studied by one or 
more of the above techniques, table I summarizes what we 
were able to learn about these systems by the light scatter- 
ing methods. For the sake of brevity, I will not discuss all of 
them in the same detail but rather highlight some selected 
aspects of some of them. 



""" 1 »MLI B * / 




FIGURE 20. Data acquisition and analysis system of 
quasielastic light scattering in time domain 
(Malvern autocorrelator) and in frequency domain 
(PDP 8/e) and those of elastic light scattering. 

a. Amphoteric Latex Particles 

Monodisperse latex particles have been utilized in a wide 
variety of fields including immunochemical assays [47, 48] 
and biomembrane studies [49] in addition to their more con- 
ventional applications as the markers and calibration stan- 
dards in microscopy and light scattering and as the model 
colloids [50]. Our interest coincides with the last instance 
where an amphoteric latex system with the well-defined 
number of charges can be invoked as a model for globular 
proteins and biomembrane vesicles relative to the ioniza- 
tion behavior of their surface groups. Recently, Homola and 
James [51] were able to prepare an amphoteric latex system 
without any added surfactants which now meets the require- 
ment of a well-defined charge number. The latex particles 
consist of three monomers, styrene (S), methacrylic acid 
(MA), and N,N-diethylaminoethyl methacrylate (DEAM), 
emulsion polymerized with persulfate at 70 °C. 

We were mainly interested in whether there exists any 
size change at different pH [52]. One of the transmission 
electron micrographs taken of the samples is shown in 
figure 21 which makes it evident that the sample has a 
relatively homogeneous size distribution. Results of the con- 
ductometric and potentiometric titrations, as shown in 
figure 22, give clear evidence that the latex particles are in- 
deed amphoteric in nature and the titration valences at the 
extreme pH of 3 and 1 1 are fairly symmetric at about 6 X 10 5 
electronic charges per particle. This was deduced from the 
titration results and the size determination, which was per- 
formed by elastic light scattering according to the method 
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TABLE 1. List of the four systems studied hy tipht .»cij«ennp irrkinqut'* 



SYSTEM 



Amphoteric Lalt-K Panicles 



TECHNIQUES 



STRl'CTl 'RE at PROPERTIES 



VARIABLES 



ELS 
QLS 
EPLS 



Surface Layer Dilation; 
electrostatic interaction* 



I H 



Photoreceptor 

Relinal Disc Membranes 



Nemaiir Liquid Crystals 



ELS 
QLS 
EPLS 



Optical Birefringence 
NMH Lin.' Splittings 
FDS 



L osmotic deformation 

2. membrani lateral stiffness 

3. inn permeability 

+. surface ion binding 



1. photochemical stales 

2. chemical potential gradient 

3. ion concentration 



1. Ork'iuaiion Di«iril>uli->ii 
Function 

2, order fluctuation dynamics 



1. min-nemalogeu etim posit inn 

2. temperature 



Linear Maeromnlecules 



r'I.» 



internal normal modes 



ELS: elastic light scattering— total scattered intensity profile 
QLS: quasielaitic light scattering— angular dependent spectral shape 
EPLS: electrophoretie |j K ht faltering-electrical field depended QLS 
FDS: forward depolarized sea tic ring 




FiWRF. 2i Transmission r j 



*lr«, micropapfc f amptmteriekle* particle. uJOStm. 
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ml of HCI (0.9871 Ml 



FIGURE 22. An example of potentiometric and conductometric titration. For the conductometry, the 
left intersection point of extrapolated lines is taken as the starting point of protonation of amino 
groups and the right intersection point as the termination point of protonation of carboxylate ions. 



outlined earlier. In figures 23 through 25, we show the scat- 
tering profiles, plotted according to eq (9), at nominal pH of 
3, 7 and 11, all at the same ionic strength of 1 mM. We dis- 
play in figure 26 all the radius data so obtained at different 
pH from 3 to 11 at the same ionic strength. It is evident that 
there is a symmetric size change with respect to pH and the 
point of minimum radius at around pH 7 coincides with the 
isoelectric point determined by titration. The Stokes radii 
obtained by quasielastic light scattering agree with those in 
figures 23 and 25, 125 ±3 nm and 123 ±3 nm respectively at 
pH 3 and 11, while that of 118 ±2 nm at pH 7 does not 
agree with the radius in figure 24, 1 12 ±2 nm. Collecting all 
these results, we propose a model as depicted schematically 
in figure 27. It consists of hydrophobic core mainly con- 
stituted of styrene monomer and of hydrophilic shell made 
predominantly of ionic comonomers, MA and DEAM. The 
observed size change is then attributed to the chain expan- 
sion in the shell layer due to electrostatic repulsive interac- 
tions while the core remains relatively intact with respect to 
pH changes in the suspending medium. What remains 
uncertain however is the difference between the Stokes 
radius and that determined by elastic light scattering at pH 
7. Whether the difference could be attributed to the 
thickening of the hydration layer at the isoelectric point 
must be examined by another technique such as NMR. 

Before closing I must remark that one could raise the 
issue of whether our scheme of size determination is indeed 
probing the outer radius of swollen latex particles as con- 
trasted to some ill-defined average of the inner and outer 
radii. After performing a set of simulation studies with con- 
centric spheres having different segment densities in the 
shell volume relative to that in the core, we are convinced 
that our method is likely to underestimate the outer radius 
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FIGURE 23. Light scattering profile of amphoteric latex at pH 2.97 where 
the polarization condition was unpolarized/unpolarized. Solid curves 
represent an isotropic solid sphere with the indicated radius in the 
Rayleigh-Gans-Debye scattering. 
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FIGURE 24. Same as figure 23 at pH 6.75. 
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FIGURE 25. Same as figure 23 at pH 11.01. 



FIGURE 26. The particle radius versus pH. The dashed curve represents 
the degree of ionization and the solid curve a model of simple polyelec- 
trolyte effect for the shell layer. 
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FIGURE 27. Schematic representation of pH induced swelling and deswelling of latex particles. 



in the swollen state. Hence the latex particle dilation could 
easily be larger than what we report here. Unambiguous 
confirmation of the concentric sphere model is not possible 
with light scatterings alone. It might be possible to employ 
small angle neutron scattering [53] with perdeuterated core 
(with fully deuterated styrene monomer) and an appropriate 
DjO/HjO mixture to contrast match the core and solvent 
such that the scattering due to the shell volume can be ac- 
centuated, whereby we could indeed measure the shell 
thickness dilation with pH. 

b. Photoreceptor Disk Membrane Vesicles 

The vertebrate visual process involves complex sequences 
of events in transducing photochemical energy to electrical 
energy [54-56]. A vast literature exists concerning the role 
of the visual pigment membrane in this process [57, 58]. Our 
goal has been to isolate disk membranes from vertebrate 
rod outer segment (ROS) as intact as possible [59] and focus 
on their static and dynamic structure relative to the photo- 
receptor function. I shall now discuss what we have learned 



about these membranes by isolating the disk membranes 
from ROS, swelling them into vesicles in hypotonic media 
and examining them in dilute suspensions by quasielastic 
and elastic light scatterings. 

A typical homodyne power spectrum is shown in figure 28 
and the spectral halfwidths obtained at different scattering 
angles are displayed in figure 29. From these, we deduced 
the Stokes radius of 0.51 ±0.05 pm. If the vesicles were 
spherical in shape, then modelling them as spherical shells 
was quite reasonable because the bilayer thickness [60] of 
about 75 nm was negligibly small compared to 500 nm for 
the Stokes radius. Analogous to eq (9), one obtains from Eqs 
(5) through (8) in the limit of p = q =0, 



x?P(d)=sm 2 x 



(23) 



for spherical shell. The results of elastic light scattering 
plotted according to eq (23) are shown in figure 30. The 
equilibrium radius calculated from the Bragg condition of 
QR = md2, where the maxima are given by rc=odd integers 
and the minima by rc=even integers, is 0.48 ±0.06 /mi. In 
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FIGURE 28. An example of homodyne power spectrum of the scattered 
light from photoreceptor disk membrane vesicles at 50° scattering angle. 
Each solid curve in the upper figure represents a single Lorentzian profile 
with a halfwidth Av,, of 16 Hz. In the lower figure, the normalized resid- 
uals of fitting to the single Lorentzian is plotted against frequency v. 
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FlCURE 29. Avw versus k 2 /n. Three independently prepared samples are 
distinguished by different symbols. 



view of the agreement between the results of quasielastic 
and elastic light scattering, we take the spherical shell 
model for the vesicles to be valid in this particular suspend- 
ing medium [61, 62]. 




sin(0/2) 

FIGURE 30A. The scattered intensity modulation profiles of two indepen- 
dently prepared vesicle suspensions of ROS membranes are displayed by 
plots of / vv (0) sin 2 (0/2) against sin (0/2). The two sets are vertically shifted 
to exhibit the reproducibility of the extrema positions. The solid curves 
x*P(x) or (QR.Y<P(6)> is plotted against sin (0/2). The parameters used for 
the theoretical profiles are R = 0.48 ^m for monodisperse shells (A = °°) 
and R w = 0.48 \im for polydisperse shells of the indicated distribution 
index h. Zimm-Schulz distribution function is used for the number fraction 
of shells with radius R as 



^■ (/?) = ir(^r)" , ' R * exp[ " 



(h + \)R 



where R. is the weight average radius. 



FIGURE 30B. Comparison of the intensity profiles of solid sphere (solid) 
and sperical shell (dashed) is given by two theoretical curves; the experi- 
mental data for polystyrene latex sphere (0.45 fim radius) are represented 
by the circles. Here, the plot is 1(0) versus sin (0/2) unlike those in (30a). 
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Upon changing the osmotic pressure of the suspending 
medium by an impermeable second solute, e.g., sucrose, we 
observe that the vesicles deform by deswelling, which is 
caused by the chemical potential gradient of the principal 
permeable component, namely H 2 0. From the osmotic 
deformation behavior of the vesicles we were able to 
estimate the lateral compressive modulus of the membrane 
bilayer [63]. This turned out to be around 3 X 10 3 Pa. 

Since the membrane vesicles respond to the chemical 
potential gradient of water across the bilayer, their defor- 
mation behavior can be used to probe the permeability of 
ionic components including the hydrogen ion. Elastic light 
scattering studies of the vesicle shape have shown that the 
hydrogen ion is completely permeable within the pH range 
of 6-8 [64]. This is shown in figure 31 where the scattering 
profile of the spherical shell shape of the vesicles is hardly 
affected. 

The binding of Ca* 2 on vesicles was then studied by elec- 
trophoretic light scattering [65]. The Doppler shift spectra 
all at 7° scattering angle and 20 °C at a constant ionic 
strength of 1 mM are displayed in figure 32. The corre- 
sponding electrophoretic mobility profile at different Ca +2 
concentrations is shown in figure 33 where the solid curve is 
drawn with a two-binding sites model. It can be represented 
by 

r/C=K l C(n l -r) + n&(±^j£-) (24) *S 



where r is the number average bound Ca* 2 per vesicle, hj 
and K x are respectively the number of high affinity, a sec- 
ond order cooperative binding sites and the corresponding 
binding constant, and n 2 and K 2 are the other set of con- 
stants for low affinity, a first order non-cooperative binding. 
By replotting the data in figure 32, a Scatchard plot shown 
in figure 34 results, from which we estimate n x = 
(1.4±0.1)XlO* and K x = (7±2)Xl0 10 M" 2 while n 2 and K 2 
are subject to a good deal of uncertainty. 

c. Binary Nematic Solutions 

Here, our interest was focused on examining how the 
phase behavior [66] and dynamic twist modes [67] were af- 
fected by mixing of non-nematogens (biphenyl and benzene) 
to a thermotropic nematic liquid crystal, methoxy benzyli- 
dene butyl-aniline (MBBA). The scattering technique is the 
FDS method where the director of the nematic system is 
oriented parallel to the polarization of incident light while 
that of scattered light is perpendicular to the director axis 
[68]. Two examples of the FDS power spectra are displayed 
in Fig 35. With use of the twist viscosity of Gahwiller [69], 
we were able to determine the twist elastic constant of pure 
MBBA as a function of the reduced temperature. This is 
shown in figure 36. Incorporating the concentration de- 



pendence of a non-nematogen to that of temperature, we 
could deduce the mapping of the twist diffusivity coefficient 
over the entire binary nematic phase region. In figure 37, 
such a mapping is displayed where the non-nematogen is 
biphenyl [68]. Fearing that this might have been unique to 
the biphenyl-MBBA system, benzene-MBBA nematic solu- 
tion was also examined [70]. The result, as shown in figure 
38, indicates that such a behavior is not restricted to the 
biphenyl-MBBA system. 
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FIGURE 31. Light scattering intensity profiles of the vesicles at different 
pH to show that hydrogen ion is completely permeable to the membranes 
within the pH range studied. Profiles give no indication of the osmotic 
deformation which is clear evidence for shape invariability. 
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FIGURE 32. Three heterodyne spectra of electrophoretic light scattering, 
all at 7° scattering angle and 10V peak-to-peak (27 V/cm) at three different 
Ca* 2 concentrations at 1 mM ionic strength solution 




FIGURE 34. A Scatched plot where the number average degree of binding 
per unit ligand concentration r/C is plotted against the degree of binding 
r. Error bars apply to every data point although only two are shown. Solid 
curve is drawn according to the model discussed in the next while dashed 
line is drawn by ignoring lower five points in r. The intercept with the 
abscissa of the dashed line gives the apparent first order binding site 
number n, and its slope the apparent first order dissociation constant K,' 1 . 
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FIGURE 33. The electrophoretic mobility, corrected to 20 °C in water, is 
plotted against Ca* 2 concentration in log scale. Error bars apply to every 
data point although only three are shown. The solid curve is drawn ac- 
cording to the model discussed in the text. 
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FlCURE 35. Heterodyne power spectra of pure MBBA at 36.90 °C, Av,^ = 
30.1±0.1 Hz (upper) and at 44.95 °C,Av 1A = 6.9±0.6 Hz (lower). 
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FlCURE 36. The elastic constants of pure MBBA versus the reduced tem- 
perature t(T). The bend (a — 3) and splay(ff=l) elastic constants are from 
the work of Haller and the twist constant (a = 2) is ours with use of the twist 
viscosity of Gahwiller. 




FlCURE 37A. Composition reduced twist diffusivity coefficient Ku/y,$(X) 
versus the reduced temperature t{T, X) = T— T S x{X)/T^(X), where T ml (X) 
= (0)(1-250,Y) as observed from the phase diagram. The filled circles are 
from temperature scan of pure MBBA, the open circles are from 
temperature scan of 0.081 mole fraction biphenyl solution and the crosses 
( + ) are from composition scan at 23 °C. 

FlCURE 37B. <t>(X), the composition dependent factor of the twist diffusivi- 
ty coefficient at constant reduced temperature versus mole fraction of 
biphenyl. 
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FIGURE 38. Same as figure 37 for the MBBA-benzene system. 



d. Intramolecular Chain Dynamics of Isotactic Polystyrene 

A typical FDS power spectrum [46] of an isotactic polysty- 
rene sample (Afu> =3.5 XlO 6 , <5 2 >'^ = 1270A) in tetrahy- 
drofuran at 4.0 mg/mL is shown in Fig 39. All spectra were 
analyzed according to eq (21) with the Zimm spacing [71] of 
relaxation times up to five relaxation modes. Truncation 
beyond the 6th mode was called for due to the chosen band- 
width of 2 kHz whereby the higher modes were obscured by 
the shot-noise level. Spectral analysis was effected by a 
3-parameter fitting routine; the parameters were the spec- 
tral intensity, shot-noise level and the terminal relaxation 
time Tj. By fixing the spacing to that of Zimm type, what we 
extract from the experiment is T lt the slowest relaxation 
time of internal normal modes. Because all measurements 
were performed at finite concentrations, we had to extrapo- 
late the data to infinite dilution. This is effected by plotting 
Tj against relative viscosity t\ t of polymer solution and ex- 
trapolating to f) r ~ 1. Such a plot is shown in figure 40 from 
which Ti at infinite dilution is determined as (1.2-1.5) XlO" 4 
sec with an experimental uncertainty of ±30%. The theo- 
retical value of 7j for the polymer sample is 1.5 XlO" 4 sec. 
Within the experimental error the two are in agreement. 



7. Concluding Remarks 

It is my hope that the readers of this report would be as 
excited as I have been about the potentials of light scatter- 
ing techniques and have some appreciation for what kinds 
of problems can be tackled. I must again emphasize that 
this is intended only as a sketchy review of what my students 
and I have been able to do and not as an up-to-date review 
of the field. I have cited several excellent reviews in the text. 
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Figure 39. FDS power spectrum of isotactic polystyrene (M.W. 3.5 X 10 6 ) 
in THF at C = 4.0 mg/mL. The solid curve is for a single Lorentzian fit 
and the dashed curve is for Zimm fit with the first five relaxation times 
(Z5). 




FIGURE 40. The terminal relaxation time deduced from the Z5 fit of the 
FDS power spectra of isotactic polystyrene in THF is plotted against rela- 
tive viscosity at different concentrations (0.65 mg/mL — 4.0 mg/mL. Error 
bars define 95% confidence limit of the fit and the theoretical value of the 
Zimm model is designated as T? tt in the ordinate. 
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